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Abstract 



We present the full 0{oi$) computation of the interference effects between the 
Higgs diphoton signal and the continuum background at the LHC While the main 
contribution to the interference originates on the gg partonic subprocess, we find that 
the corrections from the qg and qq channels amount up to 35% of it. We discuss 
the effect of these new subprocesses in the shift of the diphoton invariant mass peak 
recently reported by S. Martin in Ref. pQ. 



March 2013 



*defio@df.uba.ar 

t nfidanza@df .uba. ar 

■frogerOdf.uba.ar 

§jmazzi@df.uba.ar 

^yrotstein@df.uba.ar 

"gfsborlini@df.uba.ar 



Recent discoveries of a new particle presented both by CMS [2] and ATLAS [3] collaborations 
at the Large Hadron Collider (LHC) are consistent with the Standard Model Higgs boson. Up to 
now, its mass is estimated to be around 125.3 ± 0.4(stat) ± 0.5(syst) GeV by CMS and 126.0 ± 
0.4(stat) ± 0.4(syst) GeV by ATLAS, though the experimental uncertainty could be reduced to 
100 MeV [4] after collecting enough luminosity. 

The resonance observed in the reconstruction of the diphoton invariant mass in proton-proton 
collisions at the LHC turns out to be one of the main discovery channels and, therefore, requests 
precise theoretical calculations for the corresponding cross section. 

The dominant mechanism for SM Higgs boson production at hadron colliders is gluon-gluon 
fusion |5J, through a heavy-quark (mainly top-quark) loop. The QCD radiative corrections to the 
total cross section have been computed at the next-to-leading order (NLO) in Refs. [SHE], and 
at the next-to-next-to- leading order (NNLO) accuracy in [9HTT]. NNLO results at the exclusive 
level can be found in Refs. [T2T[T5] . State of the art computations for this channel [16] include 
electroweak corrections at NLO [TTJUH] and soft gluon resummation to next-to-next-to leading 
logarithmic accuracy 

The rare decay H — > 77 is also mediated by {W and heavy quark) loops [22j[23]. Corrections 
to the corresponding branching ratio are known at NLO accuracy for both QCD [81 I24TT29] and 
electroweak [18] cases. Missing higher orders are estimated to be below 1%. 



The corresponding background for diphoton production has also been recently computed up 
to NNLO [30] , but the interference between signal and background has not yet been evaluated to 
such level of accuracy. 

The interference of the resonant process ij — > X + H{— > 77) with the continuum QCD back- 
ground ij — > X + 77 induced by quark loops can be expressed at the level of the partonic cross 
section as: 

s:;i q/ s ^2 \R 6 (.•Aij^x+H-AH^~ n A* ont ) 



s — m 2 H ) 2 + m 2 H r 2 H 
2mHTH (s - m 2 H ) 2 + m 2 H T 2 H ' (1) 



where s is the partonic invariant mass, tuh and Th are the Higgs mass and decay width respectively 

a 

As pointed out in [H2JE3], given that the contribution arising from the real part of the ampli- 
tudes is odd in s around m H , its effect on the total 77 rate is subdominant. For the gluon-gluon 
partonic subprocess, Dicus and Willenbrock [33] found that the imaginary part of the correspond- 
ing one-loop amplitude has a quark mass suppression for the relevant helicity combinations. Dixon 
and Siu [32] computed the main contribution of the interference to the cross-section, which origi- 
nates on the two-loop imaginary part of the continuum amplitude gg — > 77. Recently, Martin pQ 
showed that even though the real part hardly contributes at the cross-section level, it has a quan- 
tifiable effect on the position of the diphoton invariant mass peak, producing a shift of (9(100 MeV) 
towards a lower mass region, once the smearing effect of the detector is taken into account. 



'For a review see [201121) . 

■I- The details on the implementation of the lineshape |3T] have a very small effect on the light Higgs discussed 
in this note. We rely here on a naive Breit-Wigner prescription. 
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Figure 1: Sample of Feynman diagrams contributing to the interference 



The gg interference channel considered in pQ is not the only 0(a\) contribution that has 
to be taken into account for a full understanding of the interference term, since other partonic 
subprocesses initiated by qg and qq can contribute at the same order. While these channels are 
suppressed with respect to the gg subprocess for the Higgs signal, they dominate the 77 QCD 
background, and therefore their contribution to the interference can not be a priori neglected. At 
variance with the gg subprocess that necessarily requests at least a one-loop amplitude for the 
background, the contribution from the remaining channels arises from tree-level amplitudes, and 
can therefore only contribute to the real part of the interference in Eq. ([T]) [§. 

In this note we present the results obtained for the remaining qg and qq channels, finalizing a 
full (lowest order) calculation of the interference between Higgs diphoton decay amplitude 

and the corresponding continuum background. We concentrate on the effect of the new interference 
channels on the position of the diphoton invariant mass peak. 

The amplitudes of the qg and qq initiated contributions to the interference were calculated 
using the Mathematica package FeynArts [33], and the analytical manipulation to obtain the 
final squared matrix element of the complete interference was done with the help of the package 
FeynCalc [35J. A sample of the Feynman diagrams for the qg interference channel are shown in 
Figure [U The diagrams and amplitudes of the remaining channels can be obtained by performing 
the corresponding crossings. 

It is worth noticing that, compared to the gg subprocess, there is an additional parton in the 
final state in the new channels. This parton has to be integrated out to evaluate the impact on 
the cross section, and its appearance might provide the wrong impression that the contribution 
is next-to-leading order-like. However, since signal and background amplitudes develop infrared 
singularities in different kinematical configurations, the interference is finite after phase space 
integration and behaves as a true tree-level contribution, with exactly the same power of the 
coupling constant as the one arising from the gluon-gluon interference channel. 

For a phenomenological analysis of the results we need to perform a convolution of the partonic 
cross-section with the parton density functions. We use the MSTW2008 LO set [36] (five massless 

§ Apart from a small imaginary part originated on the heavy-quark loop in the Higgs production amplitude that 
is discarded in this note since we rely on the effective ggH vertex. There is also an imaginary contribution in the 
decay process H — > 77 which was included since the full expression for the vertex was used, but was found to be 
negligible. 
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Figure 2: Diphoton invariant mass distribution for the interference terms. The solid line is the gg channel 
contribution, the dotted one the qg channel, and dashed the qq. 



flavours are considered), and the one- loop expression of the strong coupling constant, setting the 
default factorization and renormalization scales equal to the diphoton invariant mass (/ip = = 
Mjj). For the sake of simplicity, the production amplitudes are computed within the effective 
Lagrangian approach for the ggH coupling (relying on the infinite top mass limit), approximation 
known to work at the few percent level for the process of interest. The decay into two photons is 
treated exactly (using m t = 172.5 GeV, m b = 4.75 GeV, m c = 1.40 GeV [IBIH], m T = 1.776 GeV, 
m w = 80.395 GeV [37J) and we set a = 1/137. For the continuum background gg — > 77 we 
consider five massless flavours. For the Higgs boson we use m# = 125 GeV and Th = 4.2 MeV. 
For all the histograms we present in this section, an asymmetric cut is applied to the transverse 
momentum of the photons: p^^ so f^ > 40(30) GeV. Their pseudorapidity is constrained to 
|?7 7 | < 2.5. We also implement the standard isolation prescription for the photons, requesting that 
the transverse hadronic energy deposited within a cone of size R = A(f> 2 + Ar/ 2 < 0.4 around 
the photon should satisfy pr,had — 3 GeV. Furthermore, we reject all the events with i? 77 < 0.4. 
The effect of the precise definition of the isolation prescription is negligible since no final state 
photon-quark singularities appear at the level of the interference. Therefore, almost the same 
results are obtained if the smooth isolation prescription [38] is implemented. 

In Figure |2] we show the three contributions to the full signal-background interference as a 
function of the diphoton invariant mass M 77 after having implemented all the cuts mentioned 
above. The gg term (solid line) represents the dominant channel, while the qg contribution 
(dashed) is about 3 times smaller in absolute magnitude, but as we can observe, has the same 
shape but opposite sign to the gg channel. The qq contribution (dotted) is two orders of magnitude 
smaller than the gg one, and with the same sign of qg. The position of the maximum and minimum 
of the distribution are located near M 77 = M# ± r#/2, generating a shift which at this level 
remains at (9(1 MeV). 



3 




123 124 125 126 



(GeV) 

Figure 3: Diphoton invariant mass distribution for the interference terms including the smearing effects 
which simulate the detector (gaussian function of width 1.7 GeV). The solid line represents the gg channel 
contribution, the dashed line represents the qg channel, and the dotted one, the qq. 



To simulate the smearing effects introduced by the detector, we convolute the cross-section with 
a gaussian function of mass resolution width cr MR = 1.7 GeV following the procedure of Ref. [TJ. 
The corresponding results are presented in Figure [31 where we can observe that the magnitude of 
the interference is reduced, and the position of the peak (and dip) is moved as much as 2 GeV. 
The main reason for this shift is the highly antisymmetric nature of the interference, which is 
enhanced by the convolution with the gaussian function. The precise value of the displacement 
is proportional to both the width of the gaussian, and the absolute magnitude and sign of the 
interference, resulting in a shift towards larger invariant masses for the qg and qq channels. 

In order to quantify the physical effect of the interferences in the diphoton invariant mass 
spectrum, we present in Figure H] the corresponding results after adding the Higgs signal. The 
solid curve corresponds to the signal cross-section, without the interference terms, but including 
the detector smearing effects. As expected, the (signal) Higgs peak remains at M 77 = 125 GeV. 
When adding the gg interference term, we observe a shift on the position of the peak around 
90MeV towards a lower mass (dotted), compatible with Ref. pQ. If we also add the qg and qq 
contributions (dashed), the peak is shifted around 30MeV back towards a higher mass region 
because of the opposite sign of the amplitudes. Given the fact that qq and qg channels involve one 
extra particle in the final state, one might expect their contribution to be even more relevant for 
the corresponding interference in the process pp — > H(—t 77) + jet, since the usually dominant gg 
channel starts to contribute at the next order in the strong coupling constant for this observable. 

It is worth noticing that the results presented here are plain LO in QCD. Given the fact 
that very large K-factors are observed in both the signal and the background, one might expect 
a considerable increase in the interference as well. As a first attempt to try to quantify the 
theoretical uncertainties, we study the factorization and renormalization scale dependence of the 
interferences by varying both scales (simultaneously) from |M 77 to 2M 77 . The results are shown 
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Figure 4: Diphoton invariant mass distribution including the smearing effects of the detector (gaussian 
function of width 1.7 GeV). The solid line corresponds to the signal-only contribution. The dotted line 
corresponds to the distribution after adding the gg interference term, and the dashed line represents the 
complete Higgs signal plus all three interference contributions (gg, qg and qq). 

in Figure [5j were we present the contribution from each channel including the smearing effect 
with <7mr = 1.7 GeV. Very large variations are observed, reaching up to 45% at the peaks in all 
channels. 

Unfortunately, it is known that the LO scale variation for both signal and background in Higgs 
production is not enough to estimate the size of the missing higher order corrections. Therefore 
one might consider the result of this variation only as a lower bound on the contribution of higher 
order terms for the interference. 

While calculating the interference to the same accuracy of the signal and background seems 
hard to be achieved at the present time, a prescription to estimate the uncertainty on the evaluation 
of the interference and a way to combine it with more precise higher order computations for signal 
and background for gg — > ZZ was recently presented in [39] . The procedure can be easily extended 
to the case presented here by including all possible initial state channels. 

Finally, we would like to emphasize that a more realistic simulation of the detector effects 
should be performed in order to obtain reliable predictions and allow for a direct comparison with 
the experimental data. The Fortran code to compute the interferences presented in this paper can 
be obtained upon request from the authors. 
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Union under the Grant Agreement number PITN-GA-2010-264564 (LHCPhenoNet). 
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Figure 5: Scale dependence of the diphoton invariant mass distribution for interference terms. The gg 
channel contribution is shown in (a), qg in (b) and qq in (c). Factorization and renormalization scales 
are varied simultaneously as fj,p = fiR = kM 77 in the range ^ < k < 2. 
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